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We show that an appropriately defined fluctuation-dissipation theorem, connecting generalized 
susceptibilities and time correlation functions, is valid for times shorter than the nucleation time of 
the metastable state of Markovian systems satisfying detailed balance. This is done by assuming that 
such systems can be described by a superposition of the ground and first excited states of the master 
equation. We corroborate our results numerically for the metastable states of a two-dimensional 
Ising model. 
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There have been many efforts to extend the concepts 
and methods used to describe systems in equilibrium to 
systems which are not in equilibrium but are either sta- 
tionary, or evolve very slowly [3, |(j . A particular class 
of "slowly evolving" out of equilibrium systems are those 
which are in a metastable state and, due to their ubiquity, 
their characterization is of especial interest. Usually it is 
thought that the macroscopic properties of a metastable 
system can be treated as if it was in equilibrium. In par- 
ticular, even relations such as the fluctuation dissipation 
theorem (FDT) 0,013 

are generally assumed to be valid 
for systems in a metastable state. However, metastable 
systems are actually far from equilibrium and there is 
no reason to expect the validity of this theorem for such 
systems, even if their evolution is very slow. Indeed, 
the FDT does not apply to systems such as finite-range 
spin glasses, domain growth processes, structural glasses, 
among others (See Ref. Q and references therein). In this 
letter we use a dynamical approach to show why it is 
justified to apply results from equilibrium to metastable 
states for the case of Markovian stochastic dynamics and 
we derive the FDT for these systems from the microscopic 
dynamics. 

Since the phenomenology of metastable states has 
been assumed to be similar to that of equilibrium sys- 
tems, most of the efforts have focused on understand- 
ing the mechanisms by which a system decays from the 
metastable state to equilibrium by nucleation processes 
(growing of a second phase) 0, H, 0l El El ■ However, a 
theory for the description of metastable states per se is 
still lacking 

El 

This is partly because the phenomenon 
of metastability is a relative and rather complicated con- 
cept El- Penrose and Lebowitz El; made a detailed 
characterization of the principal properties observed in 
the behavior of systems in a metastable state, which can 
be summarized as follows: In a metastable state, a sys- 
tem behaves similarly to a hypothetical pure thermody- 
namic phase, although the intensive parameters have val- 
ues such that the equilibrium state would consist of a dif- 
ferent phase or a coexistence of different phases. When 



the system is isolated, the metastable state remains for a 
very long time. The response to small and slow perturba- 
tions leads to small and reversible changes in the systems. 
For large or rapid changes, the system may escape irre- 
versibly from the metastable state. Beyond qualitative 
characterizations, there is not a clear and general defini- 
tion of metastability |l4| . In this work we use a defini- 
tion of a metastable state similar to that introduced by 
Davies El El for Markovian systems satisfying detailed 
balance, in terms of the eigenvectors of the corresponding 
time independent master equation El- Using this defi- 
nition and the Kubo formalism for linear response theory 
HI HI ) we obtain a metastable fluctuation-dissipation 
theorem valid for times short compared with the nucle- 
ation time of the system. 

In the following we limit ourselves to the dynamics of 
Markovian systems with a finite (possible large) number 
of states. Those can be described by a master equation 
to which we associate an operator Lq 

P(t) = L P(t), (1) 

If the system is characterized by a set of discrete random 
variables a, then P(t) is a vector of components p(a, t) 
which correspond to the probability that the system is in 
the state specified by a at time t. When Eq. Q is solved 
by separation of variables, we obtain the time indepen- 
dent master equation 

L tpj = -Clji/Jj, (2) 

where flj and ipj are the eigenvalues and eigenvectors 
of L 0l respectively. Due to conservation of probability 
there exists a stationary solution ipo{&) associated to the 
eigenvalue £lo = 0, namely 4>q{<j) — p e {3). Here p e is 
the Boltzmann probability distribution of the system in 
the equilibrium state, which we assume to be the only 
stationary solution of the master equation. If detailed 
balance holds then Lq is self-adjoint with respect to the 



2 



following internal product 



R{S)Q{S) 



(3) 



with R and Q two arbitrary functions with finite 
norm[l8|. For R = tpo an d Q = V'jj Eq. implies 
^2gipj{^) = $o,j) i- e - VtoO?) is normalized and ipj^oi^) 
sum to zero. 

We assume now that it is possible to choose the pa- 
rameters of the system in such a way that one of the 
eigenvalues of Lq, labeled by — flj., corresponds to a de- 
cay that is much slower than the observational times. 
This is < ili *C 1 <C for j > 2, in appropriate 
units. This assumption means that we neglect the case 
of having several different metastable states. The exten- 
sion to finitely many metastable states is straightforward, 
but not the one to systems with a divergent number of 
metastable states (glasses, spin glasses, etc.). 

Let us now prepare the system in any configuration 
a'. Since the set of eigenfunctions is complete (as follows 
from the self-adjointness of Lq), we can represent the 
corresponding probability distribution as 



<W = ^o(^) + 



3=1 



(4) 



and hence, for I/O2 <C t <C one gets a nearly 

stationary state (which essentially does not vary in time 
for t < 1/fii) 



(5) 



where G(ct') = %l) 1 {a , )/ij) {a'). That the RHS of Eq. © is 
a probability distribution follows from the norm and pos- 
itivity preserving properties of exp(Lot). When G(a') <C 
1, the state given by Eq. (J5J is the equilibrium state, since 
the second term in the RHS is negligible. On the other 
hand, when G(ct') ~S> 1 the state is sharply localized in 
a zone of configurations {(?} m (hereafter, the metastable 
zone) and very small outside this zone. Then, this qua- 
sistatic probability distribution represents the metastable 
state. Notice that in this situation G is independent of a' 
since any configuration prepared within the metastable 
zone is expected, on physical grounds, to evolve into the 
same intermediate metastable state p m (a). The case 
G(a') ~ 1, leads to configurations which have compa- 
rable probabilities of evolving into either the equilibrium 
or metastable state. We assume that this set of "sad- 
dle points" is negligible compared to the sets of both, 
equilibrium and metastable configurations. 

Now consider a system which can be prepared in a 
metastable initial state described by 



As p m (<i) is negligible outside the metastable zone, we 
approximate p m (a) as 



p m (a) 



^o(a) + G^i(ct) for a e {a}, r 
for a $ {a} 1 



(7) 



with G> 1. The last equation coincides with the def- 
inition of the metastable state given by Davies 0, 0] , 
where the reader can find greater detail. 

Using the definition of G and the fact that it is constant 
within the metastable zone one gets tpi(a) = Gipo(cr) for 
a G {(?} m . On the other hand, using that p m (cr) = 
outside the metastable zone, we get tpi(<j) = —1/Gi/jo(<j) 
for a £ {<?} m Thus, the first excited state ipi(a), and 
hence, the metastable state, is specified in terms of 
the equilibrium distribution for (almost) all configura- 
tions of the system since it is locally proportional to 
the Boltzmann distribution in both, metastable and not 
metastable zones. The proportionality coefficients are 
given by G and — 1/G, respectively. 

This simple picture of metastability allows us to go 
beyond the description of the distributions characterizing 
the metastable states. In particular, we now derive a 
FDT through linear response theory for these states. We 
now consider the perturbed master equation: 



P + hPi = (L Q + he wt L x (P + hP x ) 



(8) 



where L\ is the perturbative term generated by an oscil- 
latory external field and P is the probability distribution 
in the absence of the perturbing external field, whose 
evolution is described by Eq. QJ. 

Now, if the system is initially in its metastable state 
— described by Eq. © — after some algebra, we obtain 
the following general expression for the changes in the 
probability distribution Pi(ct, t), to first order in h. 



P 1 (a,t)=h / e^'-'Ui [M*) + Gil>i (<?)] e lult 'dt' 
Jo 

+h [ e^-^Lx Ghj)i{&) e iuJt ' ( e - nit ' - l)dt' (9) 



This expression is exact for all times. Since L\ , tpo an d ipi 
are independent of t' the integrations are trivial. Expand- 
ing Pi (tr, t) in the basis of eigenvectors of Lq, the equilib- 
rium case is recovered for t — > 00. For times t <C l/f2i, 
the second integral vanishes, and the first integral yields 
the total change of the probability distribution starting 
from the metastable initial condition. 

Introduce for any p the following notation 



{B(t)) p = J2B(a)p(a,t). 



(10) 



p m (a) = MJ) + GMt); 



G> 1. 



(6) 



We can calculate the changes of the average value of any 
physical quantity B(a) as (B(t))p 1 where Pi is taken 
from Eq. . 
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By taking the corresponding Laplace-Fourier trans- 
form of {B(t))p i; we can define [l9| the metastable sus- 
ceptibility of the system as 



(s + Slj) 



six 



(s + Slj)(s + Ox — ilu) 



(11) 



Then, in linear approximation, the metastable suscep- 
tibility consists of two terms. The first term is similar to 
the equilibrium case but this time the initial condition is 
the probability distribution of the metastable state. The 
second is a memory term corresponding to a convolution 
with the external field. 

We now define ipj(h) as the eigenfunctions of the oper- 
ator Lq + hh\. To first order in h we have the following 
relation connecting Lx?pi(0) with the derivatives oitpx(h) 
with respect to h: 



ii^(O) = -Lo 



dh 



■fii(O) 



h=0 



dh 



h=0 



~dh 



h=0 



^(0) 
(12) 



We substitute Eq. (|12|l in the scalar products of 
Eq. (fTT|l and use the appropriate proportionality between 
ipx ((?) and (<?) ■ We then split the sum over a in a sum 
over the stable zone and one over the metastable zone 
(See Eq. EJ, noting that if the system has a metastable 
state as defined above, then the higher excited states 
satisfy T,ff m 4>j(^m) = E<?0{<? m } ~ °> for 3 ^ 2 - 

This must be the case as J2se{a } B(<t, t) — 1 f° r times 
shorter than the nucleation time 

We now define E(S) as the energy of the system ap- 
pearing in the Boltzmann distribution and introduce, for 
each probability distribution p(a,t), the following dy- 
namical correlation 



B(t) 



P a 5' 

(13) 

where p(a' , t'\a, t) is the conditional probability that the 
configuration a' occurs at time t' given that it was in a 
at time t. After several pages of algebra one then gets 
for the metastable susceptibility 



X (s) = PC 



Sli& ( B(t) — (E 



kT 



hxSl x )J + Six j dre lw{t - T) A (^B(t)^ + o 



(Si!) 



(14) 



where f3 = 1/jfeT and A(A) = {A) pm - (A) r . 

To obtain the metastable fluctuation-dissipation theo- 
rem it is enough to show that the second and third terms 
in Eq. I|14|l vanish as the coexistence curve is approached 
(first order correction in Six). Since B(t) and B(t) re- 
main bounded, the terms related to them are negligible 
because Sl\ -C 1. Indeed, the principal correction is given 
by the term of the order of d In Six I dh. Since Six is the nu- 
cleation rate|ld|. it is given roughly by exp(— @W) where 
W is the nucleation barrier. Usually W is roughly i?^ -1 , 
where R c is the critical droplet radius and diverges alge- 
braically as the supersaturation ho goes to zero. Here the 
total strength of the external field is given by the fixed 
initial field ho plus the perturbation h. (ho in the Ising 
model is the external magnetic field). From this follows 
that dSlx/dh also diverges algebraically in h, whereas Six 
goes to zero as a stretched exponential. Thus we have 
established the central result of this work, a fluctuation 
dissipation theorem for the metastable states. 

In order to check the FDT for metastable states, we 
will show that the susceptibility obtained by perturbing 
the metastable state of a two dimensional Ising model 
with a magnetic field of fixed frequency, agrees with that 
obtained by taking the Fourier-Laplace transform of the 
correlations of the fluctuations in the metastable state, 



given by the first term in the r.h.s. of Eq. I|14|) . Our 
programs where proved checking the well known FDT in 
equilibrium (See Fig. We obtained the autocorrela- 
tion of the magnetization by a Monte Carlo simulation 
for a two dimensional Ising model evolving by Glauber 
dynamics [lij on a square lattice with periodic boundary 
conditions. The set of external parameters (temperature 
T, coupling J and static magnetic field ho) was chosen 
to give long-lived metastable states (~ 10 4 Monte Carlo 
time steps per spin) when starting with all spins opposite 
to the magnetic field. As suggested in Ref. |2fJ,|2lJ, we 
used T = gT c J, where T c is the critical temperature and 
J is the coupling between nearest neighbor spins. The 
external magnetic field ho was chosen by trial. For all 
calculations we allowed the system to evolve until a long- 
lived state opposite to the field was reached. Then we 
computed the average of the fluctuations of the magneti- 
zation with respect to the metastable equilibrium value, 
over a set of 100 realizations. 

In Fig. we show the magnetic susceptibilities as a 
function of the frequency. The empty squares correspond 
to the classical magnetic susceptibility for the equilib- 
rium state, obtained by the definition of linear response 
theory. The crosses correspond to the same quantity 
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frequency 



FIG. 1: Magnetic susceptibility for a two dimensional Ising 
model with 10 4 spins. The temperature, T = 1.52, external 
magnetic field ho = 0.1 and coupling J = 1. As usual, the 
real part is the positive quantity and the imaginary part is the 
negative. The crosses correspond to the metastable computa- 
tion by the perturbative method and the empty squares to the 
equilibrium one. The solid and dashed lines were obtained by 
transforming the metastable and equilibrium autocorrelations 
of the magnetization. The inset shows the autocorrelation of 
the magnetization. The solid and dashed lines correspond to 
the metastable and equilibrium correlations respectively. 

for the metastable state. Both quantities were aver- 
aged over 10 realizations. The dashed lines correspond 
to the equilibrium magnetic susceptibility computed us- 
ing Laplace- Fourier transformation of the corresponding 
autocorrelations (per spin) of the magnetization given 
in the inset (dashed line). Finally the solid lines corre- 
spond to the magnetic susceptibility obtained by trans- 
forming the appropriate autocorrelation function in the 
metastable state (continuous line in the inset). We find 
the same good level of agreement for the equilibrium and 
metastable cases. It is interesting to observe that the be- 
havior of the correlations are very different in both cases. 

In summary, by using a formal definition of a 
metastable state for the case of Markovian systems with 
a finite number of states we have shown that the FDT in- 
deed holds for times shorter than the nucleation time. We 
also evaluated the size of the leading corrections. Since 
many systems have Markovian dynamics on sufficiently 
long time scales, this result has quite a broad range of ap- 
plicability. A crucial hypothesis was the existence of one 
single low-lying excited state of the operator Lq. This 
means that nucleation is the slowest physical process, a 
condition often satisfied in practice. Detailed calcula- 



tions and numerical simulations will be given elsewhere 
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